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The pairing properties of nuclear systems are a sensitive probe of the effective 
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nucleon-nucleon interactions. We compare the Sq pairing gaps in nuclear and neu- 
tron matter derived from the phenomenological Gogny interaction and a renormal- 
ization group motivated low-momentum V\ ow k interaction extracted from realistic 
interactions. We find that the pairing gaps predicted by these interactions are in an 
excellent agreement in a wide range of sub-nuclear densities. The close agreement 
between the predictions of the effective forces remains intact in the case where the 
single particle spectra in neutron and nuclear matter are renormalized by nuclear 
£h ■ interactions. 
> 

X 

The effective nuclear forces are the key ingredients in the nuclear structure 
studies of finite nuclei within the Hartree-Fock-Bogoliubov and related mean 
field theories. Well-known examples are the zero-range Skyrme [1] and finite- 
range Gogny [2,3] forces that have extensively been used in the large scale 
numerical calculations of finite nuclei over several decades. The effective forces 
are commonly adjusted to the bulk properties of nuclear systems after a mean- 
field variational minimization of the ground state energy of a collection of 
nuclei within the density dependent Hartree-Fock or Hartree-Fock-Bogoliubov 
theories. 

A distinctive feature of the finite range Gogny forces, which will be discussed 
below, is that these were adjusted by minimizing the Hartree-Fock-Bogoliubov 
energy functional that is a function of the pairing fields [2,3]. Thus, the pair- 
ing correlations in the system, seen experimentally in the odd-even staggering 
effects caused by the pairing in the isospin T = 1 states, are encapsulated in 
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the effective force. Since, in general, the effective forces are subject to simple 
parameterizations and have the advantage of reducing the numerical cost of 
the extensive nuclear structure calculations, it remains an important task to 
scrutinize the reliability of effective forces in different contexts, in particu- 
lar their relation to the interactions derived from the underlying microscopic 
theories. 

During the recent years much effort went into formulations of nuclear inter- 
actions in terms of effective field theories [4,5,6,7,8,9,10,11]. The main goal of 
these theories is the separation of the long-range component of the nuclear 
forces, which is dominated by the pion exchange and is well under control, 
from the intermediate and short range components, which are dominated by 
correlated pion and heavy meson exchanges that are poorly known. 

A line of approach, developed by the Stony-Brook group, applies the renormal- 
ization group arguments to the Lipmann-Schwinger (LS) equation to eliminate 
the high- momentum modes of a phase-shift equivalent potential Vnn which 
serves as a driving term in the LS equation [12,13,14,15,16,17,18]. Note that 
low momentum nucleon-nucleon (NN) interactions were also derived by the 
Bochum-Jiilich group by applying a method of unitary transformations to full 
(un-truncated) meson exchange interactions derived from chiral Lagrangians 
[10,11]. The Stony-Brook low-momentum NN interaction V\ ow k is obtained by 
integrating out the high momentum components of Vnn beyond a scale A. 
The following LS equations for the scattering amplitudes with driving terms 
Vnn an d V\ ow k are considered: 

f-OO I 

T(k',k,k 2 ) = V NN (k',k)+ q 2 dqV N N(k',q)-r- 2 ^— — T(q, k, k 2 ), (1) 

Jo Ar — q z + i\3 + 



Tiow k(p',P,p 2 ) = Mow k(p',p)+ / q 2 dqV iow k{p\q)— , n+ Ti ow k (q,p,p 2 ). 

Jo p A — q z + iu + 

(2) 

Note that the intermediate state momentum q is integrated from to oo 
and to A in the first and second equation respectively. The equivalence of 
the T-matrices derived from the LS equation above is required: T(p',p,p 2 ) = 
T\ ow k(p',P,P 2 )', (p',p) < A. As described in Ref. [14], the Vi ow k interactions 
are then derived by the Andreozzi-Lee-Suzuki method [19]. 

The V\ ow k so derived reproduces the empirical deuteron binding energy, NN 
scattering phase shifts up to E lab = 2h 2 A 2 /M. Experiments give us informa- 
tion about phase shifts only up to E\ ah ~ 350 MeV. Thus an appropriate choice 
for A is ~ 2 fm _1 . Beyond this momentum, Vnn is model dependent and lacks 
physical ground. An interesting feature of the resulting effective interaction 
^iow k is that it is largely independent of the underlying microscopic force that 
was used as a driving term in the LS equation [18]. This property clearly re- 
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fleets the fact that the low-momentum part of the microscopic interactions is 
well constrained by the experimental data, i.e., they equally well reproduce 
the binding energy of the deuteron and are phase-shift equivalent [14,18]. The 
initial applications of the Vj ow k interaction are in the shell model calculations 
[14,16,17] and in various treatments of the equation of state and the pairing 
in nuclear systems [15,20]. Clearly, the renormalization group based decima- 
tion procedure, adopted to derive the effective V\ ow k interaction, bridges the 
gap between the effective and microscopic interactions in a controlled manner, 
thus it provides a good starting point for mean-field calculations. 

The purpose of this Letter is to compare the predictions of the Gogny and 
Viow k interactions for the pairing in nuclear systems in the 1 Sq interaction 
channel. We will be concerned with infinite, zero-temperature matter, param- 
eterized in terms of the Fermi-momentum fop and will mainly focus on effects 
of the force on pairing. We shall briefly comment on the renormalization of 
the single particle spectrum in the mean-field approximation, but leave aside 
the issues of the vertex and propagator renormalizations beyond the mean 
field. Some aspects of the pairing, which are complementary to this study, 
are explored in Refs. [15,20] and [21,22] using the Vj ow k and Gogny interac- 
tions, respectively. Our work, in part, is motivated by the observation that 
the Gogny interactions predict pairing properties that are surprisingly close 
to those derived from the bare realistic (phase shift equivalent) nuclear inter- 
actions [21]. 

The effective Gogny interactions are of the generic form 



i.e. they contain two Gaussian terms which reflect the finite-range of the in- 
teraction, and a contact density-dependent term responsible for short range 
correlations. Here p is the density and P a , P T are the spin and isospin exchange 
operators. For completeness we reproduce the values of the parameters accord- 
ing to the D1S parameterization in Table 1 [3]. (Compared to the original Dl 
parameterization [2] the D1S force gives a lower surface tension and, at the 
same time, a smaller even-odd staggering which is closer to the experimental 
values.) 

The V\ ow k interaction we shall employ is based on the Nijmegen 93 potential 
with a momentum cut-off A = 2.5 fm _1 , which corresponds to a distance 
scale A -1 ~ 0.4 fm. The choice of the underlying realistic potential is not im- 
portant, as the effective interactions derived from various realistic interactions 
are practically identical (c.f. [18]). 



V( ri - r 2 ) = t (l + x P a )S(^i - r 2 )p d 
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+ E i W ™ + B m P a - H m P T - M m P a P T ] exp 



m=l 
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The gap in the quasiparticle spectrum of infinite nuclear systems is governed 
by the Bardeen-Cooper-Schrieffer (BCS) integral equation for the gap function 



A(k) = -\ f A dk'k'*V(k,k>) A(k,) (4) 

where V(k, k') is the momentum space effective pairing interaction, Sk is the 
quasiparticle spectrum, ep is the Fermi-energy, and A is momentum space 
cut-off. The high-momentum behaviour of the pairing interaction, for the case 
of Gogny force, is controlled by the finite range of the Gaussians; hence we 
can take safely the limit A — > oo in Eq. (4). The V\ ow k interactions have a 
sharp cut-off at high-momenta, and it is natural to identify A with the cut- 
off in the interaction. Note that when transformed in the momentum space 
the Gogny interaction depends on the momentum transfer in the process, and 
hence on the angle between the relative incoming and outgoing momenta of 
the particles. It is then suitable to average the matrix elements over the angle 
between the vectors k and k'. The angle averaged 1 So-wave pairing interaction 
can be written as [21] 



V(k, k') 



y/ftkk 



- C« ex P 



sinh 



fi m kk' 



(5) 



m=l 



where C m = fj, m (W m — B m — H m +M m ) . Note that only the density independent 
part of the Gogny interaction contributes to the pairing in the 1 S'o-channel. 
The angle averaged pairing interaction thus contains two Gaussian describing 
long range attraction and a short range repulsion (the terms oc C\ and C2, 
respectively) . 



Next we need to specify the single particle spectrum, in the gap equation 
(4). We shall employ two approximations. First, the single particle spectrum in 
a non-interacting limit will be used to understand the correlations between the 
pairing gap and the pairing force. Second, the single particle spectrum will be 
renormalized within the mean-field approximation. We shall use the Hartree- 
Fock single particle spectra for the Gogny interaction, which are derived below. 
For the V\ ow k interaction we shall employ the Brueckner-Hartree-Fock scheme 
(see Ref. [20] and references therein.) It should be kept in mind that when the 
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w m 


B m 


H m 


M m 


1 


0.7 


-1720.3 


1300.0 


-1813.53 


1397.6 


2 


1.2 


103.639 


-163.483 


162.812 


-223.934 



Table 1 

The parameters of the D1S Gogny interaction [3]. The ranges of the interaction jx m 
are in fm while the remainder coefficients are in MeV. The values of the parameters 
of the contact term are to = 1390 MeV fm 4 , xq = 1, and d = 1/3. 
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effective interactions depend on energy, e. g. when one is dealing with time- 
retarded interactions, the wave-function renormalization differs from unity 
and tends to counter-act the reduction of the mass caused by the momentum 
dependent self-energies. Thus, the net effect of medium renormalization of 
particle mass could be an overestimate. 

To define the single-particle spectrum, it is useful to start with the expression 
for the ground state energy (at zero temperature) in the mean-field approxi- 
mation 



h 2 k 2 1 

E = 2 ~2^ Ui + 2 Sfa' I V I ^ - 3 i )' n * n j> ( 6 ) 

i ij 

where the first term is the kinetic energy, the second term is the potential 
energy of the mean-field interaction; the indices refer to the nucleonic states, 
are their occupation probabilities. In the case where the interaction is density 
dependent the single particle potential is given by the functional derivative of 
the second term (Ei nt ) in Eq. (6), i.e. Ui = SE int /Srii, or explicitely 



1 

Uri = Ys&\ v \ V ~ i^ n 3 +i)E0 ,J l q-| 3 l - l j) n 3 n h (7) 

3 ^ L jl ° Pt 

where r(= n,p) is the isospin index, Q is the volume. Equation (7) should be 
solved self-consistently with the normalization condition for the total density 
Ei rii = p. 

The matrix elements defining the single-particle potential are evaluated in the 
mean-field approximation using plane waves for the nucleon states. We find 



pd+ i 



U T (k) = |(2 + x )(2 + d)p d+l - |(1 + 2x ) [2x T + d(x\ + xl) 
+ \ E i( 2W ™ + B m ) - (2H m + M m ) x T ] 

1 m=l 

+ EE F -( k " k ') \( H m + 2M m ) [n p (k') + n n (k')} - (W m + 2B m ) n T (k') 

m=l k' 



(8) 



where x T = p T /p are the relative concentrations of the neutrons and pro- 
tons, n T (k) = 9(kp T — k) are their occupation probabilities, and we defined a 
momentum dependent form-factor as 

F m (k - k') = 7r 3 / 2 /4 l e^l k - k 'l 2 / 4 . (9) 
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The first two terms in Eq. (8) correspond to the contributions of the contact in- 
teraction, where the terms proportional to d originate from the rearrangement 
interaction. The third and the fourth terms are the direct and the exchange 
contributions of the finite range part of interaction. In the zero-temperature 
limit of interest the phase space integrals over the Fock term in Eq. (8) can 
be done analytically, 



U* ock (k) = -- E [W m + 2(B m - M m ) - H n 



m=l 



I — 



exp -^(k + k Fq y - exp -^(k - k 



Fq) 



+Erf 



■ik + k Fq ) 



-Erf 



(k - k Fq ) 



(10) 



and the net single particle energy U T = JJ^ artrcc + U^ ock , where {/"Hartree corre _ 
sponds to the first three terms in Eq. (8), becomes an analytical function of 
the Fermi-momenta of neutrons (k Fn ) an d protons (k Fp ). The single particle 
spectrum is then defined as e T (k) = k 2 /2m + U T (k). Although the equations 
above are specified for arbitrary isospin asymmetry, we shall further concen- 
trate on the two special (symmetric nuclear matter) and x p = 
(pure neutron matter). 

To understand the correlations between the effective pairing interactions and 
the pairing gap it is useful to fix the argument of the gap function on the 
left-hand-side of Eq. (4) at the Fermi-momentum, i.e., 



A(k F ) = -I [ A dk'k' 2 V(k F ,k') 
2 jo 



A(k') 



(e k ,-e F y + A(k>) 2 



;n) 



The kernel of the gap equation is a product of the momentum space matrix 
element and V(k F ,k) which, as we shall see below, is a smooth function of 
momentum k and the anomalous propagator (the remainder multiplier in the 
kernel) which is bell-shaped with the maximum 1/2 for = e F . Since the main 
contribution to the integral comes from the vicinity of the Fermi-surface the 
differences in the pairing gaps reflect the differences in the effective forces in 
the vicinity of the Fermi- surface. This observation motivates an approximation 
where we replace V(k F ,k) by V(k F , k F ). Such an approximation permits to 
solve Eq. (11) analytically and the result is the well know BCS weak-coupling 
formula. Table 2 lists the values of the gaps calculated numerically (without 
approximations) using free single-particle spectrum along with the diagonal 
elements of the effective interactions. A clear correlation is seen between the 
ratios of the gaps derived from Gogny and V\ ow k interactions and correspond- 
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A Gogny 


Ayiow k 


R A 


^Gogny 


Mow k 


R v 


0.1 


0.05 


0.10 


0.552 


-39.5 


-49.1 


0.80 


0.2 


0.35 


0.46 


0.758 


-38.4 


-46.5 


0.83 


0.3 


0.83 


0.98 


0.852 


-36.6 


-42.9 


0.85 


0.4 


1.42 


1.57 


0.907 


-34.3 


-38.9 


0.88 


0.5 


2.03 


2.15 


0.946 


-31.6 


-34.8 


0.91 


0.6 


2.60 


2.66 


0.978 


-28.6 


-30.8 


0.93 


0.7 


3.09 


3.04 


1.006 


-25.4 


-26.9 


0.94 


0.8 


3.38 


3.26 


1.037 


-22.2 


-23.2 


0.96 


0.9 


3.52 


3.27 


1.076 


-19.2 


-19.7 


0.97 


1.0 


3.47 


3.06 


1.131 


-16.3 


-16.5 


0.99 


1.1 


3.22 


2.64 


1.220 


-13.7 


-13.5 


1.01 


1.2 


2.81 


2.02 


1.391 


-11.3 


-10.7 


1.06 


1.3 


2.26 


1.28 


1.766 


-9.3 


-8.3 


1.12 


1.4 


1.66 


0.56 


2.964 


-7.6 


-5.9 


1.28 



Table 2 

The pairing gaps for different Fermi- moment a Hf computed with the Gogny (second 
column) and the V\ ow k interactions (third column). The fourth column shows the 
ratio, i?A, of the gap values given in the second and third columns to make the 
variations with the effective interaction explicite. The fifth and sixth columns show 
the diagonal elements of the interactions Vcogny^F, &f) and Vl ow k(^F> &f) and the 
last column - their ratio, Ry. 

ing diagonal matrix elements (these ratios are related in a highly non-linear 
manner). It is remarkable that the predictions of the Gogny and V\ ow k inter- 
actions differ by no more than 10% in a wide range of subnuclear densities 
(0.4 < kp < 1.0 fm _1 ) which are relevant for the description of pairing phe- 
nomena in finite nuclei. It is also seen that the pairing gaps are a sensitive 
probe of the interactions (as can be expected from the BCS weak coupling 
result): a 20% deviation in the forces results in a factor of two deviation in the 
pairing gaps. Finally, it is worthwhile to note that, there is no such correlation 
for the realistic bare interactions, for the diagonal elements of, e.g., the Paris 
or the Reid interactions are repulsive. 

The correlations between the momentum dependent pairing gaps A (A;) calcu- 
lated numerically (without approximations) with the momentum space matrix 
elements V(kp, k), which would correspond to a pairing gap derived from Eq. 
(11), can be observed in Fig. 1. Quite generally, the momentum dependence 
of the pairing gap reflects the momentum dependence of the pairing poten- 
tial. This relation becomes explicit if one approximates the pairing force by a 
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Fig. 1. The pairing gaps in the So channel and the corresponding pairing potentials 
V(kF,k) as functions of the momentum k for several fixed Fermi-momenta fcp. 
The black and grey (green) lines refer to the V\ ow k and the Gogny interaction 
respectively. The solid, dashed, dashed-dotted and dotted lines correspond to the 
values of kp equal 0.2, 0.6, 1.0 and 1.4 fm . 

separable potential. Writing, schematically, V(k, k') = g(k)g(k') and inserting 
this form in the gap equation, one finds that the solutions are of the form 
A(fc) = Cg(k) where C is a constant. 

The V\ ow k interaction is attractive below the cut-off scale A for the Fermi- 
momenta of interest and vanishes above it; therefore there is always an as- 
sociated non-zero solution to the gap equation which likewise vanishes above 
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the cut-off scale. The high momentum tail of the Gogny interaction (k > A) 
is slightly repulsive and the pairing gap changes its sign above this scale and 
vanishes at much higher momenta; clearly, these high-momentum components 
do not affect the values of the gaps at their Fermi- surface, since the kernel of 
the gap equation is sharply peaked at the Fermi-momentum. As seen in Fig. 1 
the differences between the predictions of the effective forces for the gap func- 
tions are closely correlated with their differences in the vicinity of respective 
Fermi-momenta. 




Fig. 2. The pairing gaps in the So channel as functions of Fermi-momentum 
(density) of the matter for two effective interactions. The heavy and light lines refer 
to the Vi ow k and the Gogny interaction respectively. The solid lines correspond to 
the non-interacting single-particle spectrum, the dashed and dashed-dotted lines - 
to renormalized single particle spectrum in symmetric nuclear matter and neutron 
matter, respectively. The single particle spectra are computed in the Hartree-Fock 
theory for the Gogny interaction and the Brueckner-Hartree-Fock theory for the 
V| ow k interaction. 

Fig. 2 shows the density dependences of the pairing gaps in neutron and nu- 
clear matter (kp refers to the Fermi-momentum of nucleons and neutrons 
respectively) and the effects of single particle renormalization. The single par- 
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ticle spectra for the Gogny interactions were computed in the Hartree-Fock 
theory, as described above. For the V\ ow k interactions we used the Brueckner- 
Hartree-Fock theory with continuous choice of the single-particle spectrum. If 
one adopts free single particle spectra for nucleons, the Gogny interaction pre- 
dicts pairing gaps that are systematically larger than those predicted by the 
Viow k interaction for momenta larger than about 1 fm _1 and slightly smaller 
for momenta below this value; this picture is consistent with the behaviour of 
the potentials (see Fig. 1). The renormalization of the single particle spectra 
reduces the density of states at the Fermi-surface and, hence, the magnitude 
of the pairing gap. Independent of the chosen interaction this reduction is 
larger in the symmetrical nuclear matter than in the neutron matter, since 
in the former case the single particle spectra are steeper at the Fermi-surface 
mainly due to the tensor channel 3 Si — 3 Di interaction. For neutron matter, 
where the renormalization of single particle spectrum is mild, there is a close 
agreement between the gaps computed with the two different interactions for 
&f < 1 fm" 1 and the deviations at larger densities are not dramatic. For the 
symmetrical nuclear matter the deviations are larger, indicating that the dif- 
ferences in the single particle spectra are more important for the evaluation 
of the gap than the differences in the residual pairing interaction. 

Extrapolations of the results above to finite nuclei depend on the extent the 
high-density region contributes to the average pairing gap. If such a contri- 
bution is important, the pairing gaps predicted by the Vj ow k interaction must 
be systematically smaller that those predicted by the Gogny interaction. A 
reduction of the pairing gaps by a factor of two compared to the experimen- 
tal values was observed with realistic Argonne interaction in Sn isotopes [23], 
which suggests that this might also be the case for the V\ ow k interaction. 

To conclude, we observed a close agreement in the predictions of the pair- 
ing properties of nuclear systems by two effective interactions which have 
largely different origins - the Gogny phenomenological interaction, with pa- 
rameters fitted to reproduce nuclear properties in Hartree-Fock-Bogoliubov 
calculations, and the V\ ow k interaction, which is a renormalization group mo- 
tivated low-momentum reduction of the realistic interactions. We made explicit 
the correlations between the values of the pairing gaps and the diagonal and 
"half on the Fermi-surface" matrix elements of the effective interactions. 
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